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Abstract 

Loeb and Waxman have argued that high energy neutrinos from the decay of pi- 
ons produced in interactions of cosmic rays with interstellar gas in starburst galaxies 
would be produced with a large enough flux to be observable. Their model is reex- 
amined here and we obtain an upper limit to the diffuse neutrino flux from starburst 
galaxies. The upper limit obtained here is a factor of ~5 lower than the flux which 
they predict. Our predicted neutrino flux would be below the atmospheric neutrino 
foreground flux at energies below ~ 300 TeV and therefore would be unobserv- 
able. Compared with predicted fluxes from other extragalactic high energy neutrino 
sources, starburst neutrinos with ~ PeV energies would have a flux considerably 
below that predicted for AGN models. 

We also estimate an upper limit for the diffuse GeV 7-ray flux from starbust 
galaxies to be O(10 -2 ) of the observed 7-ray background, much less than the 
component from unresolved blazars and more than an order of magnitude below 
the estimate of Thompson et al. 
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1 Introduction 

Interactions of cosmic-ray nuclei with interstellar gas nuclei in our galaxy 
produce secondary pions. The neutral pions then decay to produce most of 
the galactic 7-rays above 0.1 GeV which have been observed [1]; the decay of 
the charged pions produces galactic cosmic-ray neutrinos [2]. 

It was pointed out 30 years ago that the distribution of high energy 7-rays in 
our Galaxy is related to the distribution of molecular clouds and very young 
hot high-mass stars in OB associations which are short-lived and explode 
into supernovae [3,4]. This association between supernovae which are likely to 
produce cosmic rays and dense regions of molecular gas led to the scenario 
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where interactions between the gas and cosmic rays then produce the 7- 
rays via the decay of the 7r° mesons produced in these interactions. A natural 
implication then would be that galaxies which were undergoing a phase of 
extremely active star formation would be likely sources of high energy 7-rays. 
Such sources are known as "starburst galaxies". Detailed model calculations 
have been made to predict 7-ray fluxes above 100 MeV for the starburst 
galaxies Arp 220 [5] and NGC 253 [6]. The predicted flux for Arp 220 is close to 
the sensitivity limit for the GLAST large area telescope whereas the predicted 
flux for NGC 253 is an order of magnitude higher than that sensitivity limit 1 
and close to the upper limit obtained from the EGRET data [7] . 

Loeb and Waxman (LW) have suggested that such hadronic processes in star- 
bust galaxies, involving cosmic ray interactions with interstellar gas followed 
by the decay of tt^'s, can produce, in toto, a large enough background of 
diffuse high energy neutrinos to be observable [8] with a very large neutrino 
detector such as Icecube [9]. LW then argue that radio observations of star- 
burst galaxies imply a lower limit on the cumulative extragalactic neutrino flux 
from starburst galaxies which is within the sensitivity range of Icecube. This 
argument is examined here and the opposite conclusion is obtained. Indeed, 
we derive an upper limit to the cumulative high energy neutrino flux from 
starburst galaxies. The diffuse flux of GeV 7-rays from the same processes 
is found to be 0(1O~ 2 ) of the observed 7-ray background, much less than 
the component from unresolved blazars and more than an order of magnitude 
below the estimate of Thompson et al. [10]. 



2 Radio Emission and the Neutrino Flux from Starbust Galaxies 

LW start with the observed synchrotron emission from starburst galaxies 
which is produced by relativistic electrons in these sources [11]. They then 
make the assumptions that (1) the presence of relativistic electrons in these 
sources implies the presence of relativistic protons, (2) the protons lose essen- 
tially all of their energy to pion production, and (3), a lower limit to the energy 
loss rate of the protons can then be obtained from the synchrotron radio flux 
by assuming that all of the electrons (and positrons) which are radiating are 
from pion decay. 

Assumption (1) is a reasonable one which is supported by observations of cos- 
mic rays in our own Milky Way galaxy. Assumption (3), viz., the "lower limit" 
assumption depends on assumption (2). However, assumptions (2) and (3) can 
be questioned because (a) the synchrotron radiating electrons may be largely 
accelerated primaries rather than secondaries related to pion production and 

1 For information on GLAST see http://www-glast.stanford.edu 
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decay as is the case in our own Galaxy, and (b) the conditions in starburst 
galaxies are significantly different from those in our own galaxy. In particular, 
starburst galaxies exhibit strong "superwinds" [14]. Such winds have signifi- 
cant dynamical effects and may disrupt magnetic fields and drive protons out 
of these galaxies before they can lose all of their energy by interacting with 
interstellar gas nuclei to produce pions. In contrast, assumption (2) of LW 
assumes full trapping of relativistic nuclei in the disks of starburst galaxies to 
the point where they only lose energy in hadronic interactions with gas atoms. 
This is in stark contrast to the situation in our own galaxy. (See footnote 2.) 

These two caveats call into serious question the argument that the radio data 
can provide a true lower limit on the cumulative diffuse flux of neutrinos from 
starburst galaxies. Let us, however, ignore them and consider that assumptions 
(l)-(3) are reasonable for obtaining an analytic upper limit for such a flux. Let 
us then accept the other estimates which lead to the ratio of injected power 
of protons to electrons at a fixed particle energy, i] p / e ~ 6 and a neutrino 
luminosity which is then related to the local radio luminosity density by 

El® v {E v = lGeV) ~ (ct H /4n)([Af(dL f /dV)} f=1AGHz (1) 

where tn is the age of the universe and ( = 3 is an evolution factor which 
takes account of the fact that starburst galaxies were more numerous in the 
past [8]. LW take the local energy production rate per unit volume at a fre- 
quency f = 1.4 GHz to be ~ 10 285 W Mpc~ 3 . Let us reexamine this value for 
f(dL f /dV)] f=lAGHz . 

The local 1.4 GHz energy production rate has been derived by LW by making 
use of an important connection bewteen radio emission and far infrared (FIR) 
emission in galaxies given in the paper of Yun, Reddy and Condon (YRC) [11]. 
That paper uses the data on IRAS galaxies to derive the local infrared lumi- 
nosity density at 60 /zm to be 2.6 x 10 7 L Mpc" 3 . This total power density is 
then used by LW to obtain the 1.4 GHz power density via a strong empirical 
correlation between the FIR and 1.4 GHz luminosity densities. 

The key difference between the result to be derived here and that obtained 
by LW is in chosing how to interpret the paper of YRC. YRC state that 
less than 10% of the local FIR luminosity density is contributed by luminous 
IR galaxies with L FIR > 10 11 M Q ; this is the component which includes the 
starburst galaxies. (Figure 11 of YRC yields an estimate of ~ 6%.) 

If we take the local contribution from starburst galaxies at 60/im to be 2.6 x 
10 6 L Q (the 10% upper limit found by YRC), using the relation (4) of YRC 
as shown in their Fig. 5(a) one finds that the component of the local radio 
luminosity density related to the starburst galaxies is at most $i.4gh z — 10 18,4 
W Hz _1 Mpc~ 3 which, when multiplied by 1.4 GHz, gives ®!agh z < 10 27 5 W 
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Mpc 3 . This value is an order of magnitude lower than the flux obtained by 
LW [8]. 2 

This is not the whole story because there is a higher relative fraction of the 
energy input from the higher relative number of starburst galaxies at higher 
redshifts. To estimate this effect, we assume that the FIR background is pro- 
portional to the integrated star formation activity rate. The fraction of the 
FIR background, k(Az) contributed by galaxies in different redshift ranges, 
Az, is obtained from Ref. [12]. Then we multiply k(Az) by the fraction of the 
FIR background contributed by starburst galaxies in different redshift ranges, 
£(Az), to estimate the mean fraction of the total FIR background contributed 
by starburst galaxies. Estimates for k and £ are shown in Table I. 

Table 1: Relative contributions to the v Starburst Galaxy Flux (see text). 



Redshift Range (Az) 


k(Az) [12] 


ttAz) 


Reference for £ 


to 0.2 


10% 


< 10% 


[11] 


0.2 to 1.2 


68% 


~ 13% 


[12] 


>1.2 


22% 


~ 60% 


[13] 



Using the results from Table 1, we estimate that 23% of the observed FIR 
background integrated over redshift is from starburst galaxies. 



3 Observability of High Energy Neutrinos from Starburst Galaxies 

The upper limit on the radio flux from starburst galaxies obtained above can 
be used to obtain an upper limit on the neutrino flux from starburst galaxies 
by using equation (1) as derived by LW. One then finds that the neutrino 
background energy flux from starburst galaxies would be at most ~ 2 x 10~ 8 
GeV cm -2 s _1 sr _1 Such a flux would be undetectable above the atmospheric 

2 One might also ask why not consider more "normal" galaxies with lower FIR 
luminosities and add them all in to estimate a higher neutrino flux? We note that 
in the case of our own galaxy, cosmic rays lose only a small fraction of their energy 
in pion producing interactions before escaping the disk, contrary to assumption (2). 
Also, cosmic rays with energies above 1 PeV (the relevant range for producing 100 
TeV neutrinos) have a differential particle spectrum cx E~ 3 2 , much steeper than 
the E~ 2 spectrum assumed by LW. This spectrum probably reflects both a steeper 
composite source spectrum and a shorter confinement time than those at lower 
energies. 
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background neutrino flux, even if equation (1) is assumed to be valid when 
extrapolated to 300 TeV and even granting all of the assumptions made by 
LW. 



Table 2: Neutrino Energy Fluxes (GeV cm 2 s 1 sr 1 ) 



v Source 


£ 2 $(10TeV) 


£ 2 $(100TeV) 


£ 2 $(lPeV) 


Reference 


Atm: AMANDA-II 


2 x 1(T 6 


7 x 1(T 8 


< 3 x 10~ 9 


[15] 


Atm (Vertical) 


7 x 1(T 7 


~2x 10~ 9 




[16] 


AMANDA-II Diff.Lim. 


9 x 1(T 8 


9 x i(r 8 


9 x 1(T 8 


[17] 


Starburst Galaxies 


< 2 x 1(T 8 


< 2 x 1(T 8 


< 2 x 10~ 8 


This paper 


AGN Cores 


5 x 1(T 10 


io- 8 


io- 7 


[18] 


AGN 


3 x 10~ 9 


3 x 10~ 8 


2 x 10~ 7 


[19] 


GRB 


5 x 1(T 10 


3 x 1(T 9 


3 x 10~ 9 


[20] 


Icecube Sensitivity 




4 x 10~ 9 


4 x 10~ 9 


[21] 



Table 2 shows a comparison of the upper limit on the flux from starburst 
galaxies given here with the atmospheric neutrino flux and with approximate 
model predictions of neutrino fluxes 7-ray bursts (GRB) and active galactic 
nuclei (AGN) along with detector array sensitivities. It can be seen from this 
table that at 100 TeV none of the extragalactic sources proposed will dominate 
over the atmospheric foreground. The table shows that the present preliminary 
upper limit on the diffuse neutrino energy flux below 1 PeV from AMANDA-II 
is ~ 8.8 x 10~ 8 GeV cm -2 s _1 sr -1 in the 10 TeV to 1 PeV energy range [17]. 
The full Icecube detector array is expected to push down to a sensitivity of 
~ 10~ 9 GeV cm" 2 s" 1 sr" 1 in the energy range 100 TeV < E v < 100 PeV 
after several years of observation. 3 Under the extreme assumption that the 
primary cosmic ray spectra in all starburst galaxies are as hard as E~ 2 up to 
energies 0(10 PeV), PeV neutrinos from starburst galaxies may be detectable 
just above the projected sensitivity of Icecube. 4 However, as can be seen from 
Table 2, above 1 PeV the AGN models predict fluxes which will be significantly 
larger than the the atmospheric foreground (expected to be < 3 x 10~ 9 ) GeV 
cm -2 s _1 sr -1 at 1 PeV), as well as the fluxes from 7-ray bursts (GRB) and 
starburst galaxies. 



3 F. Halzen, private communication 

4 Even if we make a second extreme assumption that 100% of the IR galaxies at 
redshifts greater than 1.2 are starburst galaxies, we would still predict a neutrino 
flux < 3 x 10~ 8 GeV cm" 2 s _1 sr" 1 . 
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4 Observability of Diffuse 7-rays from Starburst Galaxies 

In a follow-up paper to LW, Thompson et al. [10] have estimated the con- 
tribution of 7r°-decay 7-rays from starburst galaxies to the observed 7-ray 
background in the GeV energy range. Using an E~ 2 primary spectrum they 
get estimates of E 2 Q(E) fluxes for both 7-rays and neutrinos of 3 x 1CT 7 GeV 
cm -2 s _1 sr _1 . This is a factor of 3 larger than the neutrino flux of LW because 
they estimate that all of the electrons in these galaxies are from 7r =1= decay and 
they make the further hypothesis that a significant fraction the synchrotron 
radiating electrons which emit at 1.4 GHz lose energy by processes other than 
synchotron radiation (viz., bremsstrahlung and ionization in dense gas clouds.) 
Should the diffuse differential neutrino particle spectrum continue oc E~ 2 up 
to 10 TeV and above, the flux predicted in Ref. [10] would be more than a 
factor of 3 above the AMANDA-II limit, although no such claim is made in 
Ref. [10]. 

Our estimated 7-ray flux for the same E~ 2 primary spectrum is ~ 2 x 10~ 8 
GeV cm -2 s^ 1 sr^ 1 for 7-ray energies less than ~ 10 GeV. This is O(10~ 2 ) 
of the observed flux of ~ 1.4 x 10~ 6 GeV cm -2 s^ 1 sr" 1 determined by the 
EGRET group [22]. Above ~ 10 GeV the background spectrum will steepen 
owing to absorption from pair production interactions with the extragalactic 
ultraviolet background radiation [23]. It should be noted in this context that 
Stecker and Salamon have shown that the bulk of the observed background 
can be produced by unresolved blazars [23]. Thus, the contribution to the 
diffuse 7-ray background from starburst galaxies should be unobservable. In 
this context, we note that almost all of the observed extragalactic GeV 7-ray 
sources are blazars; no starburst galaxies have been observed. 
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